Center of mass integral in canonical general relativity 

D. Baskaran, 1 '* S. R. Lau, 2,t and A. N. Petrov 1 -* 

1 Sternberg Astronomical Institute, Moscow State University, 

Universitet Prospect 13, Moscow 119899, RUSSIA 
2 Applied Mathematics Group, Department of Mathematics, 
University of North Carolina, Chapel Hill, NC 27599-3250 USA 
(Dated: 12 September 2003) 

^ ! Abstract 

For a two-surface B tending to an infinite-radius round sphere at spatial infinity, we consider the 
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Brown- York boundary integral Hb belonging to the energy sector of the gravitational Hamiltonian. 
Assuming that the lapse function behaves as N ~ 1 in the limit, we find agreement between Hb 
and the total Arnowitt-Deser-Misner energy, an agreement first noted by Braden, Brown, Whiting, 
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and York. However, we argue that the Arnowitt-Deser-Misner mass-aspect differs from a gauge 

o ' 

invariant mass-aspect by a pure divergence on the unit sphere. We also examine the boundary 

o. 

integral Hb corresponding to the Hamiltonian generator of an asymptotic boost, in which case the 

o 

I lapse N ~ x k grows like one of the asymptotically Cartesian coordinate functions. Such a two- 



surface integral defines the fcth component of the center of mass for (the initial data belonging to) a 
Cauchy surface £ bounded by B. In the large-radius limit, we find agreement between Hb and an 
^ ; integral introduced by Beig and 6 Murchadha as an improvement upon the center-of-mass integral 

first written down by Regge and Teitelboim. Although both Hb and the Beig-6 Murchadha integral 
are naively divergent, they are in fact finite modulo the Hamiltonian constraint. Furthermore, we 
examine the relationship between Hb and a certain two-surface integral which is linear in the 
spacetime Riemann curvature tensor. Similar integrals featuring the curvature appear in works by 
Ashtekar & Hansen, Penrose, Goldberg, and Hayward. Within the canonical 3+1 formalism, we 
define gravitational energy and center-of-mass as certain moments of Riemann curvature. 
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I. INTRODUCTION 



Realization of the Lie algebra g of the Poincare group G as an algebra of differential 
operators features 

M xt = xd/dt + td/dx (1) 

as the generator of a pure boost in the x-t plane. A more remarkable realization of q arises 
in classical, canonical, and special relativistic field theory, a setting where field integrals play 
the role of group generators and Poisson bracket serves as Lie bracket. The incarnation of 
the generator (1) in this setting is [1, 2] 

M xt = - J d 3 xxT a (x, t) + tP x , (2) 

where T^ u is the Belinfante tensor, the same as the stress-energy-momentum tensor for 
simple theories, and 

P x = J d 3 xT tx (x,t) (3) 

is the x-component of the canonical field momentum. Integration in these expressions is over 
an inertial three-dimensional hyperplane £ ~ E 3 determined by fixation of the coordinate 
time t. Notice that the first term on the RHS of Eq. (2), the one built with the material energy 
density T'*(x, t), defines (minus) the x-component of the field's £ center of mass, whence 
the numerical value for M tx = —M xt obtained via evaluation of the integrals is generally not 
equal to this center-of-mass component. Via the Poisson bracket M xt generates the change 
{^(x, t), M xt } in the classical field \l/(x, t) corresponding to the infinitesimal x-t boost, and 
for points x close to the origin this change is chiefly governed by tP x in (2). So indeed this 
term must be present. 

One encounters a somewhat different situation in canonical general relativity (gr) and 
the arena of spacetimes asymptotically flat at spatial infinity. Although in general such 
spacetimes possess no group of isometries, one can nevertheless realize G as an asymp- 
totic isometry group by writing down gravitational Hamiltonians which generate the asymp- 
totic symmetries. [3, 4] We draw attention to a salient feature of the resulting description. 
Namely, the generator of an asymptotic boost — say the one corresponding to M tx above— 
has a form incorporating only the first term on the RHS of Eq. (1), or Eq. (2) for that mat- 
ter. (These remarks pertain only to boost generators in GR, as the remaining generators of 
asymptotic translations and rotations look just like they should.) Moreover, the numerical 
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value for the M tx generator in GR equals the x-component of the gravitational field's center 
of mass. Rather symbolically, the value of the said generator takes the form 



where B is a large two-surface with spherical topology, d 2 x^Ja is the proper area element 
of B, and r tt is an energy density stemming from a "boundary stress-energy-momentum 
tensor" r Q/3 . [5] 

The following heuristic argument sheds light on this salient feature. For the scenario 
of asymptotic flatness towards spatial infinity, "asymptotic" crudely means fixed-time and 
arbitrarily large radial separation from the "origin," perhaps actually a closed two-surface 
with spherical topology. In E 3 such arbitrarily large radial separation corresponds — except 
on a set of measure zero, the y-z plane — to large values of x; therefore, the asymptotic 
modifier roughly indicates that xd/dt is the dominant term in Eq. (1). The ramifications 
of these simple observations for GR stem from the following key point: the numerical value 
(determined by evaluation on a classical solution) of a canonical generator in GR is a surface 
integral at infinity, that is to say a boundary integral like the one in Eq. (4) over a two- 
surface surface B enclosing the "origin" and whose points are uniformly separated from it 
by an arbitrarily large radial distance. Therefore, the integration in Eq. (4) is quite unlike 
the integration in Eq. (2) in the following sense. The subset of B on which x is small 
is itself arbitrarily small (essentially just a "great circle" around B), suggesting that the 
expression (4) — analogous to (minus) xd/dt only — is permissible. In the end, of course, a 
test for whether one has chosen the correct boost generator is whether it serves its part in 
a consistent representation of g under the Dirac algebra determined by the Poisson bracket. 
The asymptotic boost generator we consider below has long since measured up on this 
count. [3, 4] Our simple discussion here is meant only to draw attention to the discussed 
feature, one seemingly neglected in the literature. However, we do point out that within 
the framework of the Lagrangian (rather than Hamiltonian) field formulation of GR, Ref. [6] 
considered integrals of motion at spatial infinity which were in the spirit of Eq. (2) and 
special relativistic theory, while (on classical solutions) their numerical values would be of 
the form Eq. (4). Since a Lagrangian approach does not select preferred £ hyperplanes, 
both types of integrals would thus be treated on the same footing. 

Brown and York have written down a geometric expression for the integral appearing in 




(4) 
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Eq. (4). They consider the following boundary term belonging to the "energy sector" of the 
gravitational Hamiltonian: [5, 7] 



with B and d 2 X\fa as before. Note that B is (perhaps one element of) the boundary <9E of 
a hypersurface E. N is a smearing lapse function, k is the mean curvature associated with 
the embedding of B in E, and k\ Tel is the reference mean curvature of B associated with an 
isometric embedding of B in an auxiliary Euclidean three-space E ~ E 3 . That is to say 
whether B is viewed as a surface in E or in E, it has the same two-metric a^. Hence, the 
integral (5) is essentially the difference of the total mean curvatures for the two embeddings. 
As shown in [5] , such a surface integral must be added to the smeared Hamiltonian constraint 



in order to obtain a E hypersurface functional H = + H B which is differentiate on 
the standard gravitational phase space. The H functional is differentiate provided k\ ref is 
determined solely by the B metric. Other choices for k\ rei are possible [8, 9], but except 
for one passing remark will not be considered here. In appropriate limiting scenarios it is 
known that H B agrees with the standard Arnowitt-Deser-Misner (adm) notion [10] of total 
energy for spacetimes asymptotically flat at spatial infinity, the Trautman-Bondi-Sachs 
notion [11] of total energy-momentum for spacetimes asymptotically flat at null infinity, 
and the Abbott-Deser notion [12] of conserved mass for spacetimes asymptotically anti-de 
Sitter at infinity. [9, 13, 14, 15] 

In this paper we consider the the scenario of asymptotic flatness towards spatial infinity 
(spi), in which case the E gravitational initial data (hij,Kij), spatial metric and extrinsic 
curvature, obey certain fall-off conditions specified below. If we adopt this setting and 
assume that B tends to an infinite-radius round sphere at SPI, then we may obtain physical 
characterizations of the initial data in terms of the integral (5). For instance, if the lapse 
obeys iV ~ 1 in the said limit, then H generates a pure time translation asymptotically, and 
the "on-shell" (meaning on solutions to Ti. — 0) value Eoq = lim^oo Hb of H defines the 
total E energy. Braden et al. [7] and later Hawking and Horowitz [13] have noted that such 
a definition of energy agrees with the standard ADM notion of total energy, although here we 
establish this equivalence in much more detail. We argue that the Arnowitt-Deser-Misner 




(5) 




(6) 
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mass-aspect differs from a gauge invariant mass-aspect by a pure divergence on the unit 
sphere. 

Moreover, if the lapse grows like N ~ x fc , where x k is one of the Cartesian coordi- 
nate functions near spi, then H generates an asymptotic boost, and the on-shell value 
M^ k = lim^oo Hb of H defines the kth. component of the S center of mass. [3, 4] For 
this asymptotic lapse behavior, we find agreement between M^ k and the integral introduced 
by Beig and 6 Murchadha in Ref. [4] (bom hereafter) as an improvement upon the cen- 
ter of mass integral first written down by Regge and Teitelboim in Ref. [3] (rt hereafter). 
This is our first main result. We also establish the relationship between H B and a certain 
two-surface integral which is linear in the spacetime Riemann curvature tensor. This is our 
second main result. Similar two-surface integrals featuring the curvature appear in works 
on gravitational energy-momentum by Ashtekar & Hansen [17], Penrose [18], Goldberg [19], 
and S. Hayward [20], among others. We believe our results to be of relevance for comparison 
between the standard 3+1 approach to spatial infinity and more formal treatments based 
on compactification arguments. [4, 21] The results of this work complement those given in 
the seminal Refs. [3, 4], as well as those given in Ref. [16] which mentioned that the present 
work would appear. 

II. PRELIMINARIES 

A. Fall off for metric and extrinsic curvature 

Recall the definition of a spacetime which is asymptotically flat towards spatial infinity 
given by bom. Such a spacetime possesses spatial sections on which there are so-called 
asymptotically Cartesian coordinates x k with corresponding polar coordinates (r, 9, 0). Fur- 
ther, with respect to these coordinates, the large-r perturbations Shij and 5Kij of the S 
three-metric and extrinsic curvature tensor are defined by 

ha = fij + Shij ~ fij + dij{y k )r- x + b^r' 1 ' 6 (7a) 
Kij = Otj + bKn ~ dij(v k )r~ 2 + e ijr - 2 - e , (7b) 

where < e < 1. In (7a) the flat metric is diag(l, 1, 1), v k = x k /r, the 0(1) function 
ciij is of even parity [that is to say, aij(—v k ) = aij(u k )}, and the 0(1) angular function 
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is of undetermined parity. In (7b) 0^ is the zero tensor, dij is a function of odd parity 
[that is to say, d i j( y —h /k ) = —dij(u k )], and the 0(1) angular function is of undetermined 
parity. The three-metric expansion considered by RT is of the same form, but with e — 1 in 
(7a) and (7b). Cartesian differentiation of (7) yields behavior obtained via term-by-term 
multiplication by r -1 and parity reversal on leading terms. Considering (7a), one can say 
that asymptotically Cartesian coordinates define their own Euclidean background E 3 at SPI, 
with respect to which the perturbations 5hij are defined. 
In this paper we work with the expansion 

hij = fij + Shij ~ fij + a i:j (v k )r~ l + fc^r -1-6 + c^r' 2 , (8) 

with Qj also an 0(1) angular function of undetermined parity. Hence, we are allowing for 
the possibility of a single power, say for example r -3 / 2 , lying between r~ l and r~ 2 . We could 
of course adopt a more general expansion, say 

N 

K = h + ~ u + ay ("V 1 + E h( $ r ~ l ~ £q + %^" 2 » ( 9 ) 

q=l 

with up to iV terms lying between the r -1 and r -2 orders. In this case we would have 
0<£i<£2<'''<£jv^1 and bff of undetermined parity for each q. However, adoption 
of such a more general expansion needlessly complicates the analysis. Furthermore, we 
believe that all of our main results are also valid for the expansion (9). 

Let us comment on the meaning of the asymptotic symbol ~ in this paper. In the 
asymptotic expansions we write down, if the last term written is 0(r p ) with p some integer, 
then the next unwritten term is 0(r p ~ 7 ) where < 7 < 1. We have switched from e to 7 
now in order to emphasize following. For the next unwritten 0(r~ 2 ~ 7 ) term in an equation 
like (8), the 7 need not be the same as the e in the third term on the RHS. Furthermore, in 
an equation such as (7a) the ~ indicates that the next unwritten term is 0(r~ 2 ), that is to 
say the next integral power of 1/r following r -1 ~ e . Sometimes we include the order symbol 
O and write equalities for the sake of clarity. 
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B. bom integral and asymptotic scenarios 

Beig and 6 Murchadha [4] define a center of mass integral associated with the Hamiltonian 
generator of an asymptotic boost. Namely, 

M^(N) = ~^j d 2 x^n k [Nh ij (djhi k - 8 k h l3 ) + (h l k ti j - h\tii){h l3 - f^N} , (10) 

where § is shorthand for lim^oo f B with B a level-r two-surface and n k its outward- 
pointing normal. Furthermore, here, as it should for a boost [3], the lapse behaves as 

N = p ± k x k + a~ + 0{r~ £ ) , (11) 

with the P ± k = —13k 1 - constants and a~ an 0(1) angular function of odd parity. The ex- 
pression above (10) is written in terms of asymptotically Cartesian coordinates and ordinary 
partial derivatives, but it is easy to render the expression covariant. 

Were the lapse to go only as N = 1 + 0(r~ £ ), then the BoM integral would be the ADM 
energy. Indeed, one recognizes the first term in the integrand (10) as the familiar ADM 
expression, 

E 00 = -^-l d 2 x^n k h l \d 3 h ik - d^) . (12) 

We examine two scenarios in this paper; these being the energy scenario, in which the 
lapse behaves as N ~ 1, and the center-of-mass scenario, in which the lapse behaves as in 
Eq. (11). 

We note that the integral (10) is naively divergent, although its actual finiteness is ensured 
by the particular choice (11) of lapse N, the even parity of aij built into the asymptotic 
E metric (7a), and on account of the integral Hy, of the corresponding constraint H. As 
Beig and 6 Murchadha show in Appendix C of their Ref. [4], addition of M^(N) to 
given in (6) yields a total expression which is explicitly finite. (We also analyze the integral 
ifs in our Appendix C, isolating its divergent contribution via a method different than the 
BoM one.) Therefore, the BoM integral is finite on-shell, that is to say, given the vanishing 
of the Hamiltonian constraint 7i and in turn the volume term ifg. Let us use Mg{N) to 
denote the integral corresponding to (10) before the r — > oo limit is taken. The asymptotic 
expansion 

M^{N) ~ ^M^r + ( " 1+e) M^r 1 - 6 + °M^r° (13) 
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for the integral (corresponding now to a large but finite B two-surface) particularly eluci- 
dates the issues at hand. The parity conditions built into Eqs. (7a, 11) automatically ensure 
that the constant _1 M^ in fact vanishes. However, invocation of the Hamiltonian constraint 
is required to ensure that the coefficient ( _1+£ )Mg vanishes. With rt fall-off the integral 
(10) is explicitly finite without appeal to the Hamiltonian constraint, as the r l ~ e order in 
the (13) is absent. Whatever the specified fall-off is, in the end M^(N) = °M&(N). 

C. Curvature integral 

The other aforementioned two-surface integral, written in terms of the Riemann tensor, 
we wish to consider is the following: 



In this equation dt^xK is the spacetime Riemann tensor, A is the area of B, and the square 
root factor outside of the integral is asymptotic to a single power of the coordinate radius r. 
Furthermore, the spacetime representation of the B two-metric a^ u serves here to project 
free indices into B. With the future-pointing normal of S in spacetime denoted by and 
the out-pointing normal of B in E denoted by n M , we have a^ v = g^ u + u^u v — n fJ- n u '. As 
for a physical interpretation of the integral (14), consider timelike geodesies whose tangents 
at points of the hypersurface £ are given by u^. In the absence of matter u^u^^i^xuK = 
—2u fl n x u u n K ?R. tl xvK, and the latter quantity controls the rate of change of the vector joining 
any two nearby geodesies and describing the radial displacement between them. Hence, 
the integrand in (14) dictates the radial component of the geodesic deviation vector, as 
one would expect for a "mass-aspect." [21] We show below that H B and Mr also agree 
in the r — > oo SPI limit, again for both energy and center-of-mass scenarios. The on-shell 
finiteness of both H B and M^ as r — > oo in the weaker bom fall-off setting is intimately 
tied to the fact that these are correct surface integrals to add to (6), thereby achieving a 
functionally different iable total Hamiltonian H = + H B - On the other hand, Mr in (14) 
has no particular relation to the canonical 3+1 Hamiltonian; therefore, its on-shell finiteness 
is far from obvious. See Ref. [20] for a discussion of how an integral very similar to Mr stems 
from a dual-null 2 + 2 Hamiltonian description of general relativity. 

Our comparison of H B with Mr leads to the following alternative definitions for energy 




(14) 
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and center of mass for gravitational initial data sets. For the energy we find 




(15) 



where the superscript 3 means take the coefficient of the 0(r~ 3 ) term in the radial expansion 
of the curvature component a^ v 'o- XK $t,^\ UK and here the averaging is over the unit sphere. This 
definition may be easily generalized to a "superenergy" — the charge integral corresponding 
to a general time supertranslation at SPI — by placing an angle-dependent smearing function 
in front of the mass-aspect \a^ h '<J XK ^t^\ VK \l '4, in other words using iV = 0(1) rather than N = 
1 in Eq. (15). We stress that the ADM mass-aspect (determined from the integrand lifted 
from the expression for the ADM energy) differs from this manifestly gauge-invariant mass- 
aspect by a pure divergence on the unit sphere. See the concluding remarks in Section Ffor 
further details. Arguably the ADM mass-aspect does not define a valid superenergy, since 
such a definition not would generally vanish for Minkowski spacetime. For the coordinates 
M^ k of the data set's center of mass, we find 



with a similar meaning for the superscript 4. 

III. H B AND THE BEIG 6 MURCHADHA INTEGRAL 

In this section we consider the Hb integral (5), with B a level-r surface, showing that 
its r — > oo limit coincides with the bom integral (10). To set up and compute the limit, 
we shall introduce a background metric on E via a certain embedding of B into an 
auxiliary Euclidean three-space £ ~ E 3 . All quantities associated with the background will 
be denoted by either san serif or boldface letters. In particular, we use k to represent k\ rci , 
or more precisely an approximation to k\ rcf of sufficient accuracy to compute the limit. 

A. Key identities 

Let us first lay some groundwork necessary to obtain useful identities for both k — k 
and N(k — k). Assume that we have a three-dimensional slice S, equipped with two dis- 
tinct proper Riemannian three-metrics: hij with compatible covariant derivative Di and 




(16) 
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with compatible covariant derivative Dj. Viewing as the physical metric and as the 
background metric, we set 

h lk = h lk + Ah lk , h lk = h lk + Aq lk , (17) 

where the terms Ahi k and Aq lk contain all orders of perturbations. 

Assume that a radial coordinate s foliates £ into a nested family of smooth closed two- 
surfaces. Later, s will be the radial coordinate r = (x k x k )~ l l 2 determined by asymptotically 
Cartesian coordinates x k . We use s at this point to emphasize the fact that here our 
calculations are not tied to a particular type of two-surface. Loosely, we use the letter B 
to represent both the s-foliation of S itself and a particular slice (or leaf) of the foliation 
determined by setting s equal to a constant value so- Respectively, let a a b and cr ab denote the 
metrics induced on a generic level-s slice by the and metrics. Our central assumption 
is that the metrics a a b and cr a b agree in a sense made precise below. The background 
metric introduced in (17) will not be arbitrary, rather it will be defined by an (essentially) 
isometric embedding of B into Euclidean space and as such will be the flat metric belonging 
to a Euclidean three-space £ ~ E 3 . 

Let n k = Md k s be the outward-pointing normal covector for B as a submanifold of 
the Riemannian space (S , hij). Likewise, let = Md k s be the outward-pointing normal 
covector for B as a submanifold of the Riemannian space (E , h^). The function M = 
[h^ (d i s)(djs)}~ 1 / 2 ensures that n k n k = 1, and likewise M ensures the n fc n fe = 1. Above and 
in what follows, indices on physical objects are lowered and raised with and its inverse 
h tj , whereas indices on background objects are lowered and raised with and its inverse 
h* J . Again, we let An 1 = n l — n l stand for all orders of perturbations. In S coordinates the 
induced (projection) metrics for the s foliation of our two Riemannian spaces are 

a ik = h ik _ n i n k ? a ik = h ik -n i n k . (18) 

Now our requirement that the metrics induced on B by and agree can be rewritten 

as 

Aa ij = 0, (19) 

where Aa' tk = a lk — <j tk . Actually, this condition is too strong. Below, when s is r = (x k x k ) l l 2 
we shall only demand that Aa tj = 0(r~ 2 ~ e ). However, as we shall retain Aa tj in our 
calculations, there is no harm in considering this stronger agreement for the time-being, 
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and doing so more clearly demonstrates the reason for working with the inverse two-metric. 
We stress that while working with three dimensional indices, one must enforce equality of 
the inverse two-metrics as above in order to ensure that a ab = <r ab on B (the same, of course, 
as o~ a t, and cr ab agreeing on B). Indeed, if we choose s as the first coordinate x 1 , then the 
index a on the B coordinate x a runs over 2, 3 and we have 

a lk = (T lk = 0. (20) 

That is to say, o %k and a %k have only 2,3 components (i. e. B components). In contrast, 
o\k 7^ 0, <7\k 7^ and a± k ^ <7\ k in general on B. Therefore, Oi k and cr^ need not agree on 
B. In other words, the equation with lower three dimensional indices which is analogous to 
Eq. (19) is generally not valid, although, of course, a a b = cr ab on B. 

With the assumptions spelled out in the preceding paragraph, we now collect the promised 
identities. Let N be a smearing function and 

k = -D i n i , k = -D i n i . (21) 

We then have the following identities: 

k — k— \n l h kl (D k Ahu - D t Ah kl ) + n^D, (r^An®) (22a) 
+ ln l Ah u D k Aq kl + ln t D k {An l An k ) 
— |/i fe 'An*DjA/i M — |A(T ife Djn fc , 

N(k - k) = \Nrfh kl (D k Ah u - D t Ah M ) - \ {n k h il Ah ik - n l h ik Ah tk ) D t N (22b) 
+n fc D, (N^ l An k ^) + lNn t Ah il D k Aq kl + iiVn,D fc (An* An*) 
-iNh^ArjDiAhki + \ (n l h lk An l An k - n t An l An l ) D t N 
-\NAa ik D ink + In'hifcA^DjiV, 

with the first identity following from the second upon assuming that iV is constant and 
unity. We derive these identities in Appendix A.l. Note that Acr y could be replaced in 
these identities via use of the appendix Eq. (Al). 

B. Construction of and various coordinates 

We construct a diffeomorphism between S and S ~ E 3 as follows. Take a level-s surface 
B in E, say the one determined by s = s , and embed it in E. At this point we make no 

11 



assumption that this embedding is isometric. In £ assume that B is also a level-s coordinate 
surface of value s . Label the points on this level surface in £ by their coordinate values 
x a inherited from B in E. Now extend the the coordinates off of B to a system (s, x a ) 
on £ in a region surrounding B. One way of doing this would be to construct Riemann 
normal coordinates. The construction described gives us a diffeomorphism, as (s, x a ) label 
points in both E and £. Further, this diffeomorphism identifies level-s surfaces in E with 
corresponding ones in E, providing us with the set-up in Section III. A where we can work 
with a single E equipped with two distinct proper Riemannian three-metrics. We stress that 
with this construction is a flat Euclidean metric, although it need not be the trivial metric 
diag(l, 1, 1). Let us make a few comments here meant to highlight the exceptional nature 
of the foregoing construction that occurs when s is the radius r = (xi t x k ) 1 ^ 2 stemming from 
asymptotically Cartesian coordinates. Suppose that the embedding into E of the original 
two-surface B defined by s = s were an exact isometry, which we can guarantee via rather 
mild assumptions on the Ricci scalar 1Z of B. Then one would not expect that level-s 
surfaces in E neighboring this initial surface would also be isometric to their counterparts 
in E. Indeed, were this the case, one would have a foliation of (an annulus of) flat Euclidean 
space in which an infinite number of slices were exactly isometric to slices belonging to a 
foliation of the non-trivial Riemannian space (£,/ijj). This would seem to us an overly 
restrictive situation to achieve. However, while working with the coordinate r, we shall find 
it possible to nearly achieve this situation by relaxing the requirement that the isometries 
are exact, instead assuming that they hold approximately through some appropriate order 
in the small parameter 1/r. Even subject to this relaxation, k and k\ Tei will agree to an 
accuracy sufficient to compute the r — > oo limit of (5) with N(k — k) in place of N(k — k\ ). 
Equation (17) has been viewed in the system (s,x a ) of coordinates just discussed, 

h ij (s,x a ) = Ms,<) + A/^-(s,0 . (23) 

The splitting above depends on (i) the initial choice of B two-surface through the embedding 
of B in E and on (ii) how the coordinates are extended off of B once the embedding is carried 
out. (Let us just loosely say that the splitting depends on the B embedding.) Nevertheless, 
our calculations are covariant on E, and we can go to any other arbitrary coordinates. For 
example, from the system (s,x a ) on E, we may transform to a truly Cartesian system X k 
on E. Via the constructed diffeomorphism, the system X k may also be placed on E. Since 
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by definition coordinate transformation to the system X k makes hjj the Kronecker delta (we 
denote this diagonal flat-space metric by f^), adoption of the system X k on S yields the 
splitting 

h ij (X k )=f ij + Ah ij {X k ), (24) 

which is quite similar to the type of decomposition used by RT and bom in the system x k of 
coordinates, namely Eq. (8). 

Now assume that the coordinate s is in fact the radius r stemming from asymptotically 
Cartesian coordinates x k . Then we may write Eq. (23) as 

hij{r, 9, 0) = h^(r, 9, 0) + Ah tJ (r, 9, 0) . (25) 

As mentioned, now the situation will be that all large level-r surfaces will have essentially 
the same intrinsic geometry a ab , whether induced by h; L j or by h^. We show this in the next 
subsection. Notice that under the transformation (r, 6>, 0) — > x k in Eq. (25), one does not 
recover Eq. (8), of course. Rather one obtains 

h ij (x k ) = h ij (x k )+Ah ij (x k ), (26) 

where Kj{x k ) ^ in general. Eqs. (25) and (26) are in fact the same unique splitting with 
the flat background X defined by the B embedding, only the coordinates differ. Eq. (8) 
is a metric splitting with respect to the flat space E 3 defined by asymptotically Cartesian 
coordinates, while Eq. (26) represents a different decomposition into background and per- 
turbation parts, one defined with respect to the different Euclidean space X. Later we will 
have need to consider the fall-off for Ahij(x k ) in (26) as r = (x k Xk) 1 ^ 2 — > oo which is also 
defined in terms of Shij through the as yet unknown expansion 

h tJ (x k ) = f tJ + 5h lJ (x k ), (27) 

with respect to in x k coordinates. Now, it is not evident that the fall-off for 5hij(x k ) 
coincides with Shij(x k ) in (8). However, as shown below, these fall-offs are qualitatively the 
same. 

C. Isometric embedding of a slightly deformed two sphere into a flat space 

In this subsection we solve the problem of removing a distant large-r two-sphere from 
an asymptotically flat slice £ and isometrically embedding it into a flat space X. Our 
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method of solution is a perturbative one. Solution of this problem defines the transformation 
between asymptotically Cartesian coordinates x k belonging to E and Cartesian coordinates 
X k belonging to the flat space X of the embedding. 

Dropping the lowest order term, we rewrite the line-element (8) associated with the 
general E metric in terms of the polar coordinates (r ,8 , 0) associated with x k . Then, 
fixing r = r , with r ^> 1 some large constant, we define a large two-surface -B(r ), which 
we refer to as a slightly deformed two-sphere. From the E line-element rewritten in polar 
coordinates, we may obtain the line element for our slightly deformed two-sphere, 

~ r ^1 + ^ + rl+£ Jd» +r ^ ^ + rl+£ J d6d<f> 

+ ^0(l + 2 ^ + ^)tf. (28) 

It is easy to conclude that under the parity transformation V(x k ) = —x k the spherical 
coordinates behave asV(r ,0 ,<f>) = (r , n — 9 , + n). Then, keeping in mind that the Oy in 
(8) are of even parity and the coefficients in (28) are held as independent, one sees that a 
and P are of the even parity, whereas 7 is of odd parity. 

Next, we write down the line-element for the Euclidean space E in the corresponding 
spherical polar coordinates, 

dsl = dR 2 + R 2 (dQ 2 + sin 2 6d$ 2 ) . (29) 

Our goal is to isometrically embed B(r ) defined by (28) into the flat space E with line- 
element (29). We posit the existence of a coordinate transformation with asymptotic form 



R/r „ l + m£ + Wg. (30a) 

e ^ e + f XM> + ^Ml ,30b) 

, h(e,(f>) H(e,4>) , , 

which features as yet undefined angular functions. Substitution of the transformation (30) 
into (29) yields another expression for the flat metric consistent with bom fall-off. We again 
fix r = r in the resulting expression to define a two-surface B(r ) and demand that its 
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corresponding line-element 



ds 2 
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,1+e 







+ r 2 



2 ( dg 2 dh 

7 {M + sm9 d9 
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F+ m) 

2 /<9G 



d£ 2 
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+ sm 2 — 







«9# y 



r 2 , sin 2 6* 



l + -[f + gcote + — 
r V d(p 



+ 



,1+e 



F + Gcotfl + 



<9F 



•o \ d< P 
matches the line-element (28) just given. Such matching is tantamount to solving the 
isometric embedding problem. As l/r is here a small parameter, we shall -after peeling off 
an overall r 2 factor- explicitly consider the r^ 1 and r^ 1 ' 6 orders and implicitly consider the 
r~ 2 order. Hence our solution of the isometric embedding problem will only be approximate. 



Balancing terms at order r 1 , we obtain the system 



J 89 



f + gcot9 



dh 



dg . 9 n dh 



a 



P 



7- 



(32a) 
(32b) 

(32c) 



while for the case e < 1 a similar balance of terms at order r 1 6 



yields the system 



^ dG 



F + Gcot9 + 



dH 



dG 2 dH 



A 
B 

Q. 



(33a) 
(33b) 

(33c) 



In the case e — 1 the system is at order r^ 2 and the RHS of each equation in (33) is modified 
by addition of known functions of /, g, h and their derivatives: A i— > A + A'(f,g,h), 
B ^ B + B'(f,g,h), Q Q + Q'(f,g,h). Of course, we would assume that /, g, h were 
obtained via resolution of the first system (32) before turning to solve the modified second 
system. In the case e < 1 we may also obtain a system at order r^ 2 , yet a third system and 
one of the same type. This is possible because we have included the extra terms 6^ and c^- 
in Eq. (8) and therefore use a more detailed expansion than the one bom use. Considering 
all of these systems together, including the possible third unwritten system at order r^ 2 , we 
may solve the embedding equations asymptotically to a high enough order to ensure that 
o~ab — Cab = 0(r~ £ ), in turn giving the following crucial result: Aa' lj = 0(r~ 2 ~ e ). 
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We may consider our isometric embedding problem solved (to the required accuracy) if 
there is a combined solution to the systems (32) and (33) as well as the third system we 
just mentioned. In Appendix A. 2 we examine the system (32), finding a solution assuming 
simple compatibility conditions for the coefficients in (28). The remaining systems are 
formally similar and may be examined analogously. Note that conditions of regularity on 
the coefficients a, (3, and 7 (intrinsic B geometry) are expected and should mirror conditions 
on the full metric a a b which are necessary and sufficient for existence of a (suitably unique) 
solution to the full embedding equations. If the Ricci scalar (twice the Gaussian curvature) 
of B is everywhere positive (as is the case here), the task of isometrically embedding B in 
Euclidean space is Weyl's problem, a classic embedding problem of differential geometry in 
the large. In the most robust formulation of the problem known to us, that due to Heinz 
[22], existence of such an embedding is guaranteed if the B metric coefficients are twice 
continuously differentiable. Uniqueness of the embedding (up to Euclidean motions) follows 
from the Cohn-Vosson theorem. [23] 

First performing the trivial transformation from the system (i?, 6, <3>) to the corresponding 
Cartesian coordinates X k , we now seek, subject to the assumption of BoM fall-off, the 
behavior of the transformation between X k and x k . In fact, the transformation (30) can be 
rewritten as [see (A31) in Appendix A. 2 and the text thereafter] as 

X k ~ x k + °£ k {v l ) + £ i k r- £ + 1 i k r~ 1 , (34) 

whence 

dX k /dx i ~ <f\ + A k i {i/)r- 1 + B k i r- 1 - £ + CV~ 2 • (35) 

The angular functions °£ fc (z^) and A k i(u l ) are of odd and even parity respectively, while the 
angular functions e £ fc , 1 ^ fc , B k i, and C k i are not of definite parity. Note that the A k i and B k i 
here are not the A and B appearing in the system (33). 
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D. Asymptotic expansions for Ahij, Aq^ , and An 1 

With respect to the bom metric splitting (8), we first collect the following standard and 
easily derived formulae: 

hij = fij + Shij ~ fij + dif" 1 + bijr" 1 ' 6 + Cijr~ 2 (36a) 

h ij = fj + §q ij „ fj _ a ij r -i _ Wh-- 1 - 6 - (c* J - a\a ll )r- 2 (36b) 

M= 1 + 5M ~ l + la^r" 1 + \b uv r- 1 - £ + \{c vu -aja iu + \al v )r- 2 (36c) 
n*= i/* + 571* ~z/ i + i(z/ i a w -2aV)r- 1 + i(z/ J 6 w -26V)r- 1 - £ 

+ | [^(cw - a/a fcl , + |a^) - 2c\, + 2a* fe a fej , - a\a vv \ r~ 2 (36d) 

with the equalities indicating that Shy, Sq 1 ^ , 5M, and 5n* contain all orders of perturbations. 
Moreover, on the rightmost side of the equations all indices have been raised or lowered 
with the trivial flat metric and a subscript v indicates contraction with v % (for example 
a\ = a ik u k ). 

With hij(X k ) = fij and the expansion (35), we find the sought-for expansion (27) of 
BoM-type for the background metric, 

hij ~ f { j + aij(v k )r- 1 + bij-r -1-6 + c^r" 2 , (37) 

where now all sub dominant terms are pure-gauge and given by 

3ij — *^-ij ^-ji (38a) 

hij = B tj + Bji (38b) 
Cjj *A.ki~A j + Cij + Cji . (38c) 

Note that a^- has even parity, while both and c^- are not of definite parity. From the 
expansion (37) we obtain identities for h^, h tj , M, and n l which are token identical to those 
given in Eqs. (36). Moreover, these mirror identities also have the same parity behavior 
as those found in Eqs. (36). Substitution of Eqs. (8) and (37) into Eq. (26) and these 
considerations show that 

Ah^ = 8hij — 5hjj (39a) 
Aq ij = Sq ij - 5q lj (39b) 
An* = <5n*-<5n* (39c) 
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each possess precisely the same fall-off and parity properties as Shij, 5q 13 , and Sn\ although 
of course they do differ term-by-term from these. 

E. Energy integral 

We turn now to the detailed comparison between the Brown- York and BoM integrals. 
If we naively insert the identities (22) in terms of A perturbations into the Brown- York 
integral Hb, then we arrive at expressions appearing to be the ADM and BoM integrals. 
However, one cannot make such a direct comparison, as we have expressed H B in terms 
of A perturbations with respect to reference Euclidean space X, whereas the BoM integral 
is expressed in terms of 5 perturbations with respect to the Euclidean space E 3 defined 
by the asymptotically Cartesian coordinates. An elegant way to achieve the comparison 
would appeal to the technique of gauge (inner) transformations [24] based on the theory of 
continuous transformations [25] . However, we adopt a more cumbersome but straightforward 
approach. 

Let us first consider energy scenario. As demonstrated later in the text around Eq. (86) 
below, the proper area-element for a level-r two-surface B obeys 



where again dVt denotes the area-element of the unit sphere. Only the leading behavior of 
this equation is relevant for the energy scenario at hand; however, the next-to-leading order 
is relevant for the center-of-mass scenario, since a lapse N growing like r will "sample" 
this next-to-leading order. Hence we keep it here for future use. Our results from the last 
subsection further imply that 



Therefore, as far as all energy scenario limits and even center-of-mass scenario limits are 
concerned, we may view d 2 X\fa as the area-element induced either by hij or h,j. 

Since the area-element grows like r 2 and for now iV ~ 1, we see that only leading 0(r~ 2 ) 
terms in the expression (22a) for k — k will contribute to the r — ■> oo limit of Hb- It is readily 
seen that the last four terms in (22a) can not contribute to the limit. Indeed, symbolically 
D = d + r, where T Christoffel terms are 0(r~ 2 ), a fact stemming from the derived fall-off 
(37) for h^. Therefore, D k differentiation -like d k differentiation- drops fall-off by one power 



d 2 x^ = dtlr 2 [1 + \{a\ - a vv )r- x + 0{r~ x - e )\ , 



(40) 



- ^ = 0(r~ £ ) . 



(41) 
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of inverse radius (and reverses parity). This shows that the third, fourth and fifth terms are 
0(r~ 3 ). Finally, our perturbative solution to the embedding equations has shown the sixth 
and final term to be 0(r~ 3 ~ e ). 

Now considering the remaining two factors (the first two) on the RHS of (22a), we drop 
some more subdominant terms and write 

k - k ~ \tfh kl (d k Ahu - diAh kl ) + v k d x ( n ['An fe l) . (42) 

The product of the last factor and the leading contribution to the area-element (40) inte- 
grates to zero via Stokes' theorem. Whence for the energy scenario we have that 

H B~-^~J B d 2 x^n l h kl (d k Ah tl - d t Ah kl ) , (43) 

a result which yields 




after some replacements of n % and by n* and h u at the expense of introducing subdominant 
terms which are then discarded. Apart from the fact that the integration is over the finite 
surface B and not the sphere at infinity, the first integral on the RHS is the ADM energy 
belonging to E. Likewise, in the limit the second integral is the ADM energy belonging 
to the vacuum slice £ ~ E 3 . To see this, merely replace y/a with y/cr at the expense of 
introducing subdominant terms. Because we work with physically reasonable bom or rt 
fall-off, we expect that this second integral then vanishes in the limit. That it indeed does 
may also be verified via direct calculation. With Eqs. (35,37,38) we get 

K ~ fa + o, % + a, % . (45) 

Whence to leading order the integrand of the second integral in Eq. (44) is z/c^cfy ^]), and 
again the leading order of the product of this term with the area-element integrates to zero 
via Stokes' theorem. 

F. Center— of— mass integral 

Evaluation of the integral (5) for the scenario with iV behaving as in Eq. (11) is quite 
a bit more subtle. The integral is naively divergent in the r — > oo limit, so we must keep 
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track of both divergent terms (which only vanishes upon invocation of the Hamiltonian 
constraint) and finite terms. Let us turn now to the identity (22b) for N(k — k) and first 
dispatch the terms which most obviously do not contribute to the limit. First, as the last 
two terms are 0(r~ 2 ~ £ ) they do not contribute. The fourth, fifth, sixth, and seventh terms 
are each 0(r~ 2 ). Nevertheless, none of these terms contribute to the limit, since for each 
parity conditions ensure that all leading terms integrate to zero. For example, consider the 
fourth term, \Nx\ l AhuD k Aq kl . As argued just above, D fc differentiation drops fall-off and 
reverses parity, so D k Aq kl is 0(r~ 2 ) and leading odd-parity. The product rtAhu is clearly 
0(r~ x ) and also leading odd-parity. Since the lapse here is 0(r) and leading odd-parity, 
we see that the full term is leading odd-parity, whence the product of this term with the 
leading contribution to the area-element integrates to zero. 

Focus now on the third term, r\ k Di(Nn^ l An k ^), on the RHS of (22b), writing it as 

n k Di(Nn [l An k] ) = n k di(Nn [l An k] ) + n k (Nn [p An k] )T l pl . (46) 

The term involving the Christoffel symbol is 0(r~ 2 ) and with leading odd parity, essentially 
because the T k p i are themselves 0(r~ 2 ) and of leading odd parity. Therefore, this term does 
not contribute to the limit. The first term on the RHS also makes no contribution. To show 
this, first recall the discussion above and introduce the odd-parity term 

A 1 = \{v l a vv - 2a\ - v l z vv + 2a 1 ,) , (47) 

so that Art ~ A l r~ l . Next, expand Nn^An k \ the factor within the parenthesis, as follows 

Nn v An k] = /3\ v [l A k] + \ [lk] + p [lk] + 0{r'^ £ ) , (48) 

where fi L v = ft^jvi is O(r ) and and are respectively 0(r~ £ ) and 0(r _1 ). From 
the last equation, infer that 

nkd^N^An^) = n k did3\ v^A^) + + p [lk] ) + 0( r - 2 - £ ) . (49) 

With the analog of (36c) for M in n k = Mu k , the fact that fi^ v is 0(1) and of odd parity, and 
Eq. (47), we conclude that on the RHS the first and only worrisome term includes 0(r _1 ) 
and 0{r~ 2 ) pieces (both potentially contributing to the integral) which are of odd parity 
and thus integrate to zero. 
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We drop those terms in Eq. (22b) shown to make no contribution to the limit, thereby 
reaching 

Hb ~iLJ Sx ^° ^ Nnihkl ( dhAhu ~~ dtAhkl) ~ ( nhhilAhik - n W fcA M d i N ] ■ ( 5 °) 

To get this last equation we have also replaced Dj differentiation with <9j differentiation, a step 
easily seen as permissible via fall-off and parity arguments. We now write Ahij = 5hij — 5hij, 
in order to cast the last equation into the form 

Hb ~ lfL / [ NniHkl ( dhha ~ dlHkl) ~ ( nkhU5hik ~ ri hi " Sh ik) diN] 

~ iLL d2x ^ [ jvn<hW (9 * ha ~ ~ ( nfch<,<jh< * _ ■ ( si ) 

To reach this equation, we have made several swaps of h 1 - 7 for h 1 ^ and n l for n*. Such swaps 
are permissible, since they result in the introduction of 0(r~ 2 ) terms in the integrand which 
are then seen to integrate to zero via parity arguments. One recognizes the first integral as 
the bom integral M^(N) belonging to S. It is easy to infer that the second integral on the 
RHS of Eq. (51) must vanish in the r — > oo limit. First replace yfa with y/a 1 , a step which 
does not affect the limit. The resulting integral M^(iV) is the bom integral for and the 
vacuum slice S ~ E 3 ; hence we also expect this integral to vanish in the limit as h^- obeys 
physically reasonable fall-off. 

Let us quickly establish the vanishing of bom integral when evaluated on the flat metric 
hij belonging to the vacuum slice S. Now, the full integrand, the sum of four terms, is 
0(r~ x ) and of leading odd parity. Consider the product of its leading term with the area- 
element (40). It is only this leading term of the integrand which "sees" the next-to-leading 
term in the area element; that is to say, the product of the leading term of the integrand 
with the next-to-leading term of the area element potentially contributes to the integral. 
But we see that this potential contribution integrates to zero via parity arguments, whence 
in the calculations to follow we may work solely with the leading round-sphere term dVtr 2 
of the area element. 

Without affecting the r — > oo limit, the integrand for the bom integral of may be 
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replaced by the sum of the following four terms: 

terml = (^tx 1 )^ f ij djh ik (52a) 

term2 = -(0 ± ix l )v k f ii d k h ij (52b) 

term3 = - (3 1 x v { j h SKj (52c) 

term4 = ^ x v l f j 5Kj ■ (52d) 

To see why this is the case, consider, for example, the term 

NnWdjKk = Nv k f i *d j h ik + N6n k f i *d j h ik 

+ Nv k 5q ij djh ik + NSrfSj'djhik . (53) 

On the RHS the last term is 0(r~ 3 ) while the second and third terms are each 0(r~ 2 ) and 
off odd parity, whence none of these terms contribute to the limit. Finally, the first term on 
the RHS is 

Nv k pd 3 Kk = {^^ypd.Kk + a-v k fd 3 Kk + 0(r" 2 - £ ) , (54) 

where the middle term on the RHS is 0(r -2 ) and of odd parity, so on the RHS the sole 
limit-contributing term is the first one. 
Eqs. (35,37,38c) allow us to write 

h - = fij + + (E.P.T.)r- 2 + 0(r- 2 - £ ) 

= fij + 9& + d& + (E.P.T.)r- 2 + 0(r- 2 - £ ) . (55) 

Again, (E.P.T.) stands for generic terms of even parity, and the one here stems from the 
AkiA h j in Eq. (38c). When coupled with simple parity and fall-off arguments, this expression 
for hjj shows that we may replace <5hjj terms in Eqs. (52) with £%fij. Therefore, up to terms 
not contributing in the r — > oo limit, 

terml + term2 + term3 + term4 
= 2(P ± l x l )u k pd i d^ k] - >' ////''(r7,0 + d&) + 2/3V/ ij 'di& 
= 2u k d i [(P\x l )pd^ k] ]-2p ± l u i f l ^d j ^ + 2p ± l u l rd^ j 
= 2v k d i V)3 [i £ fc] - 2(3^i k] ] , (56) 

whence — to the relevant order — the full integrand integrates to zero via Stokes' theorem. 
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IV. H B AND THE CURVATURE INTEGRAL Mr 



In this section we again evaluate the main surface integral (5) in the large-sphere limit 
towards SPI, this time showing that 

HB ~~\L\f^J B d2x ^ N ^^ K ^ ( 57 ) 

as r — > oo. We have discussed the integral on the RHS in the paragraph after Eq. (14) 
where it is denoted This result holds for both the energy scenario with iV ~ 1 and the 
center-of-mass scenario with N given by Eq. (11) discussed in Section II. 
We will verify (57) by establishing the following results. First, 

k-kr { ~ 3 [yVV]r 2 , (58) 

where the superscript 3 means take the coefficient of the leading 0(r~ 3 ) term in a radial 
expansion of the curvature component ^cr^a^^^x^. The coefficient 3 [|o" m1 V Ak 9^a^k] hap- 
pens to be of even parity. Second, Eq. (40) can be written 

d 2 x ^ = r 2 dVt [1 + (E.P.T.)r- 1 + 0(r _1 - e )] , (59) 

where we use E.P.T. to stand for generic terms of even parity. Important in itself, this result 



also shows that ^A/An ~ r + 0{r~ v ). Third, 

k-k\ Tcf ~ (E.P.T.) r~ 2 + {z+e) [\^ v a XK ^ VK ] r~ 2 ~ e + (E.P.T.+ 4 [\a^a XK ^ VK ] )r- 3 . (60) 

The E.P.T. coefficient sitting before the r~ 2 is again 3 [jcr^a^^t^xvi^ , and the term 
( 3+£ ) [|o" mj/ o" Ak 3? m a^k] is n °t °f definite parity. When coupled with by now standard argu- 
ments, these three results establish Eq. (57) subject to the assumptions of either the energy 
scenario or center-of-mass scenario. 

Before turning to detailed calculations, let us lay some groundwork and fix more notations. 
We may write the £ metric in an ADM-like form 

h ij dx i dx j = M 2 dr 2 + a ab (dx a + W a dr) (dx b + W b dr) , (61) 

where M and W a are respectively the radial lapse function and radial shift vector. Take 
x<1 — (C> C) as coordinates ona5 surface. Here ( — a single complex coordinate on B — is the 
stereographic coordinate e 1 ^ cot (0/2) belonging to the asymptotically Cartesian coordinates 
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x k , and C is the complex conjugate of (. Introduce a complex null dyad m a d/dx a and 
fh a d/dx a on B, chosen so that m a rh a = 1 and m a m a = 0. In terms of the null dyad the B 
metric is written as 

Vab = m a m b + rh a m b . (62) 

One advantage of working with a complex null B dyad is that results derived for one frame 
leg e§ = m a d/dx a yield corresponding results for the other leg e~ A = rfi a d/dx a under complex 
conjugation. (See Appendix B.l for an explanation as to why we use 3 and 4 as the name 
indices here.) Another is that the components m a m b k a b and 2m a fh b k a b respectively capture 
the trace-free and trace pieces of the extrinsic curvature tensor k a b associated with the radial 
foliation of E into nested B two-surfaces. As an orthonormal co-triad on E take 

e h = Mdr (63a) 
e 3 = rh a {dx a + W a dr) = rh c d( + rh^d( + W^dr , (63b) 

where Wm = m a W a . We may obtain the one-form e 4 from e 3 via conjugation. The spatial 
triad dual to and the connection coefficients uj^ determined by this co-triad are listed 
respectively in the appendix Eqs. (B28) and (B31). Triad indices j,k,l,--- run over the 
values h, 3, 4. 

A. Asymptotic expansions 

Recall that u k = x k /r, and with the stereographic coordinate define the complex direction 

^ = v^- 1 (l - C 2 , -i(l + C 2 ) , 2C) (64) 
discussed in Appendix B.l. The Cartesian components \i % obey the relationship 

p = (jpj + + v i v i . (65) 
With fi l and u\ we define the following O(r ) quantities (and their complex conjugates): 

duu = v l v :> aij (66a) 

a-vii, = v l H j aij (66b) 

= (SpPciij (66c) 

a w = fi l fi J aij . (66d) 
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Both a vv and are in fact of even parity, although this is not obvious for the latter 
quantity, while a vll and a m are not of definite parity. Below we discuss the behavior of a v ^ 
a^p, and a w under the parity operation. 

The expansions for the metric functions M and W m are 

M ~ 1 + \a vv r- x (67a) 
W m ~ a^r- 1 . (67b) 

To motivate similar expansions for and m^, let us first note that 

Hidx 1 = rV2P' 1 d(, (68) 

which may be demonstrated via the chain rule. Therefore, fj,^ = \iidx 1 /d( = and /x<? = 
Hidx l /d( = ry/^P -1 , although both dx l /d( and its conjugate are 0(r). These results 
determine the leading behavior in the expansions for the components of the B co-frame, 



m c ~ $-r + sj\p- l a^ (69a) 
m c - ~ v^P _1 r + sj\p- l a^ . (69b) 

We have carefully tailored the next-to-leading order behavior in the expansions (69) to 
ensure that 

m b = fi b + l (a^fi b + a^Ha) r' 1 + 0(r~ £ ) . (70) 

In this equation note that fi a = ^dx 1 /dx a is 0(r). Now construct a bc via Eq. (62), thereby 
finding 

v bc = f bc + a bc r- 1 + 0(r 1 - £ ), (71) 
where f bc = ^bfic + fib^c is the 0{r 2 ) metric on a radius-r round sphere and 

a bc = (iSfib + flUb) (rfp-c + i&H^Oij (72) 

is the 0(r 2 ) projection of into the said radius-r sphere. Note that dx l /dx b = [i l ji b + 
as can be shown by Eq. (65). We see therefore that our expansions (69) for the null co-frame 
are consistent, since contraction of Eq. (7a) upon the 0(r 2 ) projector {dx % / dx b ){dx^ / dx c ) 
also yields Eq. (71). 

Let us obtain the asymptotic expansion for the action of 5 on an s spin-weighted scalar. 
Using Eqs. (69) in the appendix Eq. (B28b) for the dyad leg e§ or via simple inference, one 
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finds the not unexpected result, 

m b = / - \ (a^ b + a^fi b ) r' 1 + 0( r - 2 - £ ) , (73) 

where, for example, \i a = f ab fib = 0(l/r). Now define the operator 5 = m a d/dx a , and 
let 5 = r\i a djdx a be the corresponding operator on the unit sphere. On scalars, these 
operators agree with the ones discussed in Appendix B.l. Contraction of Eq. (73) on d/dx b 
then gives 

5 ~ 5 Q r~ l - \ (a^A + a^5 ) r~ 2 . (74) 

Next, we insert the co-frame expansions (69) into the appendix result (B31b) for the B 
connection form uj = u^, thereby reaching 

u ~ ujor' 1 + \ (5 [a^] - S [a^] - u a^ - u Q a^) r" 2 , (75) 

where c^o = — C/ V% is the connection form on the unit sphere. Now, by definition 9 = 5 — su, 
whence we find 

5 ~ 5 r _1 - \ [a^do + a w 5 + s(S a M/z ) - s(6 a w )] r~ 2 (76) 

as the sought-for 9 expansion. To get the corresponding expansion for 6, simply complex 
conjugate (76) and then send s — > —s. 

Let us obtain expansions for the dyad components of the B extrinsic curvature tensor. 
As a connection form, the trace is given by the formula k = 2k^i = — luj^-i- So we insert 
the expansions (67), (69), and (76) into appendix formula (B31c) for obtaining 

k ~ -2r _1 + (a vv + + c* a ufi + B a^)r -2 (77) 

as the desired expansion. Similarly, starting with k^ = —^31-3 and Eq. (B31d), we get 

kmm ~ + 9oa^)r -2 (78) 

as the expansion for the trace-free part of k a b- 

Finally, let us obtain the asymptotic expansion for the scalar curvature 1Z of B. Into the 
appendix formula (B34) for TZ, we insert the expansions (74) and (75), with result 

K ~ 2r" 2 + (5 2 a^ + B 2 a w - 25 9 a^ - 2a^)r- 3 . (79) 
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Our expansions for k = 2k mm , k mm , 5, and 7Z determine the leading order behavior in 
radial expansions for some components of the £ Riemann tensor. Indeed from the Gaufi- 
Codazzi-Mainardi embedding equations, [8, 26] 

(krnfh) k mm k mm 2^ Rmmmfh (80a) 
^mm S^rani Rhmmm j (80b) 

we may infer the explicit expressions for the coefficients 3 R mmm m and 3 Rhmmm m the expan- 
sions 

Rmmmfh ~ Rmmmfh^ (81a) 
Rhmmm ~ RhmmmT ■ (81b) 



Eq. (80a) determines that 



Rmmmm 2 ^ ^ ' (^) 



whence we obtain 

3 Rmmmm = —O-uu — ^O^up, ~ §0^^ — ^0 a pp — \^Q a ^^ + ^O^O^fMp • (83) 

Likewise, Eq. (80b) determines that 

Rhmmm 9(5 ^mm ^0 ^mm i (S4) 

whence 

3 Rhmmm = \ a vn — J^Oduu + — J^O^^p, ~ J^o^up + |§oSo a ^ • (85) 

To reach the last equation, we have appealed to the commutator equation 2[8 , 9o] a ^t — a u^ 
discussed in Appendix B.l. 

The following consistency check affords some confidence in our results. First, show that 

= i(|m c f - |m f | 2 ) = i2r 2 p- 2 [l + a^r' 1 + O^ 1 ^)] , (86) 

where the i is necessary for d 2 x^/~a to be real as d( A d( is pure imaginary. Now the area 
form on a radius-r round sphere is 

ifx a fx b dx a A dx b = i2r 2 p- 2 d( A d( ; (87) 
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and with this fact and Eq. (79), we perform a direct calculation showing that 

/ d 2 x^n = 8rr + 0(r- 1 - £ ) . (88) 

JB 

That is to say, the GauB-Bonnet Theorem holds to the same level of accuracy as our ap- 
proximations thus far. We point out that the presence of the sole term in 3 TZ which is not a 
unit-sphere divergence, namely — 2a M ^, plays a crucial role in this agreement, as it cancels a 
similar such term in y/a. Note that our discussion here has also established the last section's 
Eq. (40) as well as Eq. (59) at the beginning of this section. Do note, however, that the 
coordinates used here are (C, C)> whereas in Section III the B coordinates used were (9,<t>). 
These two systems are related by a non-trivial (in fact imaginary) Jacobian. Therefore, d 2 x 
here is not the d 2 x from Section III, although of course d 2 Xyfa is the same here as there. 



B. Parity 

Define the action of the parity operator V via V(x k ) = —x k , whence it follows that 
V(v k ) = —v k . Moreover, one can show that V(() = — 1/C, and from this result and (64) 
that 

W = -(C/C)^- (89) 

We can then immediately write 

V{a vv ) = a vv , V{a Vjl ) = ((/()a ufl , (90a) 

V(a^) = , V{a^) = {CKfa m . (90b) 

Moreover, 

V(5 ) = (C/C^o , V(lo ) = -(0' 2 u , (91) 

where the first of these formulae follows since P^o) = V{y\i % djdx l ). Note that the minus 
sign difference between the formulae for V(5 ) and V{\x % ) stems from V(d/dx l ) = —d/dx\ 
With the formulae amassed so far, it is fairly easy to establish the following: 

Vi^la^p) = BqCW ( 92a ) 

V(d a ufl ) = B a UfM (92b) 

V(d a^) = (C/C) 2 §oa^ (92c) 

Vidodoa^) = 8 cVmm- ( 92d ) 
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These results and their complex conjugates in turn show that 3 7Z, 2 k, 3 Rmmmm, and 
2 k m m 2 kfhm are all of even parity. 

C. Energy integral 

Let us return to the surface integral (5) and total energy scenario. Consider the expression 
(82) derived from the Gaufi-Codazzi-Mainardi equations. A similar formula is associated 
with the isometric embedding of B in S and the Ansatz 

fc|ref ^ _ 2r -l + 2 A .|ref r -2 _ (gg) 

Namely, 

2 k\ rci = -± 3 ft. (94) 

We give a more careful derivation of this result in the next subsection. Therefore, appealing 
to Eq. (82) we get 

k-k\ rei ~- 3 i? m— r" 2 , (95) 

and, noting that 

^Rmmrnfh = O ^ O R%kjl , (96) 

we may also write this as 

k-k\ rel ~ *[\a ij a kl R ikjl ]r- 2 . (97) 

Moreover, for the E components of the spacetime Riemann tensor, one has [26] 

^ijki = Rijki + K ik Kji - K u K jk . (98) 

Since = 0(r -2 ), the terms quadratic in Kij are in fact 0(r~ 4 ) terms, whence we have 
established Eq. (58). 

D. Center of mass integral 

When considering the surface integral (5) in the center-of-mass scenario, we keep more 
powers of inverse radius in the k expansion than we kept in Eq. (77). We now write 

k ~ -1r~ l + 2 kr~ 2 + (2+c)kr- 2 - £ + 3 kr~ 3 , (99) 
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with the previously given expression for 2 k still valid. From the embedding equation (80a) 
we infer that 

o(2+e)j, _l(3+e) 7? _ (3+e) d M 00^ 

2 3i, _ f2jL \2_2ii, 2l _147-._4d /-, mh \ 

" ^mm \ n mm) n mm, n mm 2 - LL mmmm j ^ iuuu i 

and so — with the results of the parity subsection — we have 
k 2r~ l + (E.P.T.)r^ 2 

_(l(3+e) n + (3+e) R _ _ \ -2-e _ /gp-p 1 14-p , 4^> _ _ \ -3 QQl) 
\ 2 L mmmm ) 1 y±j.± . j. . ~r g ' v - 1 L mmmm J 1 ■ \ u / 

Again, E.P.T. stands for generic even parity terms. We claim that 

Af cf ~ -2r- x + (E.P.T.)r" 2 - (|( 3+£ ^)r^ £ - (E.P.T. + § 4 ft)r~ 3 , (102) 

and will establish this result below. We note in passing that the expansion (102) agrees with 
the "lightcone reference" —\fTlZ derived in Ref. [9]. Combination of these equations with 
Eq. (96) gives 

k - k\ rcf ~ (E.P.T.)r- 2 + ( 3+£ ) [±a ij a kl R ikjl ] r~ 2 ~ e + (E.P.T. + 4 [\a ij a kl R ikjl ] . (103) 

Now, since Kij ~ dij(v k )r~ 2 with dij of odd parity, the quadratic extrinsic curvature terms 
in (98) are of leading r~ 4 order and even parity. Hence, we have 

(3+£) [X^^W] = (3+£) [\cr ij cj kl R ikjl ] (104a) 
4 1\<j» v <j Xk U» Xvk ] = 4 [\a ij a kl R ikjl ] + E.P.T. , (104b) 

and this result along with (103) establishes Eq. (60). 

Let us now verify (102). We must consider the GauB-Codazzi-Mainardi equations asso- 
ciated with the isometric embedding of B in £ ~ E 3 . Using k a b = (/c| rcf ) a 6 as a shorthand, 
these are 

(k m m) ^mrrMmm 2^ = ^ (105a) 

5k mm - 5k m?fl = 0. (105b) 

We solve these equations — in the sense of asymptotic expansions — for k mfh and k mm , viewing 
the intrinsic B geometry, both 3 and 71, as fixed. We make the following Ansatze: 

k m fh ~ — f + k m fj{r + ^ + ^kjnmT + k m fj{r (106a) 
k mm ~ 2 k mm r~ 2 . (106b) 
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We remark that consistency of the Ansdtze adopted here follows from the fact that - 
as shown in Section III — the non-trivial flat background metric \j admits an expansion 
analogous to hij. Therefore, one expects all of the asymptotic expansion calculated in Section 
IV. A to carry over for hy. However, here we endeavor to follow an independent approach. 
From Eqs. (106) we first algebraically find that 

2 k mm = (107a) 
{2+£) k mm = -| (3+£) ^, (107b) 

establishing in particular the claim made earlier in Eq. (94). Clearly then — as the previous 
parity analysis shows — P( 2 kmm) = 2 k mm , which justifies the E.P.T. term next to the r~ 2 
on the RHS of Eq. (102). Also algebraically, we obtain 

o\ _ - (\ -^ 2 - 2 k 2 k ! 4 7? M081 

where — as seen from Eqs. (105b) and (107a) — the weight-two scalar 2 k mm solves the pde 

B \ mm = -^ 3 n. (109) 

Lemma. V( 2 k mm ) = (C/C) 2 ( 2 kmm), whence 2 k mm 2 k mm is of even parity [which is the remain- 
ing needed piece to verify Eq. (102)]. 

To start the proof, let us show that there exists a unique solution 2 k mm to the pde given 
in Eq. (109), where we view the RHS as a prescribed source determined by Eq. (79). As is 
well-known, for an equation of the form So/ = g with g a prescribed weight-one source, a 
unique inverse 5q 1 to the operation B exists if and only if [27] 



ain 1 Y lm g = 0, (110) 



where again § dVL denotes average over the unit sphere S 2 and the \Yi m are spin-1 spherical 
harmonics. Using Eq. (B25) from Appendix B.l, we see that for our equation the issue at 
hand is whether or not 

r dn^Y lm d 3 n (in) 



vanishes for m = —1,0,1. Simple integration by parts will not show that the integral 
vanishes, since 9 (-i^im) 7^ 0. However, for any scalar curvature IZ, there exists a weight- 
two scalar Q such that dTZ = dQ. [28] More precisely, letting ip denote the conformal factor 
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relating a a b to the line-element ds^ of the unit sphere, as in 

ds 2 = 4> 2 a ab dx a dx b , (112) 

Tod shows that [28] 

Q = 4g 2 log^-4(51og^) 2 . (113) 

Assume ip ~ r" 1 + 2 ^r" 2 and that 2 ?p is of even parity. Since the log terms are all differen- 
tiated we may write 

log^ = -logr + log [l + 2 ^r~ 1 + 0(r _1 ~ £ )] 

= -logr + 2 iljr~ 1 + 0(r- 1 - £ ) , (114) 

whence (113) in tandem with the expansion (76) for 5 leads to 

Q ~ 45 2 ( 2 ^)r- 3 . (115) 

This is a consistent result since QQ = <3R = 0(r~ 4 ), the leading 2r~ 2 term in 1Z being 
annihilated by the angular derivatives in 5. We have then the leading-order identity 5o 3 7^ = 
45oc3o( 2 V ; )> an d so — dropping a factor of 4 — the integral (111) just above becomes 

£dn_ 1 Y lm o i b i l( 2 ij) = 0. (116) 

The last equality follows by integration by parts and 8o(-i^i,m) = 0. Thus, the uniqueness 
of 8 _1 is established. Whence we have 

K mm — 4°0 °0 /C 

= -a8(V). (H7) 

Since 2 ip is an even parity scalar function, the lemma then follows by simple calculations as 
outlined in the parity subsection. 

V. CONCLUDING REMARKS 

Our investigation has shown the mass-aspect to be the following: [cf. Eqs. (15,83, 96,98)] 

\ V VAK ^W] = I ( 2a - + 2g oa^ + 28 <v + + §o<W - 23 3oa^) • (118) 
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[Upon proper averaging over the unit sphere, that is to say (47r) _1 § dfl integration, the 
mass-aspect yields the total energy] This mass-aspect may be compared with the one 
stemming from the ADM integral (12). Namely, 

\ 2 [n k h^(d 3 h tk - dkhij)] = \(2a vv + cV^ + %a v „) . (119) 

(118) and (119) differ by a pure divergence on the unit sphere. We note that the reference 
term k\ mi in Eqs. (5,93) plays a crucial role in yielding a mass-aspect (118) whose proper 
unit-sphere average agrees with the ADM energy. Indeed, not only does k\ ref remove the 
leading order divergent contribution to the "unreferenced energy" [the proper B integral of 
(87r) -1 A; alone which blows up like r], at the next order it removes a dangerous factor of 
from 2 k in Eq. (77). Indeed, the total energy (Sn)^ 1 § d£la vv is not equal to (87r) _1 § dQ 2 k 
in general. 

In stark contrast with the energy scenario, we point out that the results of Appendix C.l 
(in particular the lemmas as they pertain to coefficients of the TZ expansion) show that the 
k\ rei term in Eqs. (5,102) in fact makes no contribution at all to the r — > oo limit in the 
center-of-mass scenario! This suggests that the bom integral might be compared directly 
to the proper B integral of (87r) _1 iVA; as an alternate way of checking the correspondence 
between Hb and Mg in the r — > oo limit, and one which bypasses the issue of the reference 
term and the solution to the embedding equations altogether. (Using Section IV techniques, 
we have performed such a check. The calculation amounts to a tedious exercise in perturba- 
tion theory.) However, this alternate way requires that one somehow know in advance that 
the reference term makes no contribution to the center-of-mass limit, so the reasoning would 
seem circular. Furthermore, we note that our asymptotic solution to the embedding equa- 
tions presented in Appendix A. 2 and the analysis of Section ///justify the Ansdtze (106) of 
Section IV, ultimately showing that we have needed to use this solution and that analysis in 
verifying that the reference term does not contribute to the center of mass after all. By way 
of comparison with the last sentence of the preceding paragraph, we mention that — knowing 
the reference term makes no contribution — the kth center-of-mass Cartesian component 
is indeed (Sn)~ l § <if2 3 /cz/ fc , and we find (Sir)~ l § dVL [c vv + 2c^ + 5oc^ + ^c u ^)u k as the 
explicit expression for this component. Moreover, under changes in Cy induced by coordi- 
nate transformations on hij of the form (34,35), we find that this unit-sphere integral is an 
invariant. 
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These considerations also highlight the difference between the Brown- York and bom 
integrals. Indeed, Eqs. (58,103,104) show explicitly that even the integrand of Hb in (5) 
vanishes to a high order in 1/r for trivial initial data. Actually, of course, by definition 
k — k\ rcf is identically zero to all orders if hij is a Euclidean metric and £ is E 3 . However, 
as we have seen in Section III.F, the bom integrand can be non-zero even for trivial data, in 
which case the vanishing of the full integral relies on the integration itself. This difference 
may well be relevant when supertranslations are brought into play, and we hope that our 
extremely careful treatment of the reference term will prove useful in future investigations. 
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Note added in proof: As we have often noted, our calculations have been carried out with 
our more detailed version (8) of Beig-6 Murchadha fall-off [4]. Recently, Szabados has 
carefully analyzed the Poincare structure of asymptotically flat spacetimes and found Beig- 
6 Murchadha fall-off to be the weakest possible fall-off which ensures finiteness of the angular 
momentum and center of mass [35]. 

APPENDIX A: KEY IDENTITIES AND EMBEDDING EQUATIONS 
1. Derivation of key identities 

Let us now establish the main identities given in Eqs. (22a, 22b). We begin by collecting 
some preliminary identities needed to get the main ones. Subtracting the two equations in 
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(18), we find 

Aa ik = Aq lk - n l An k - An l n k - An l An k , (Al) 
which upon contraction with n fc yields 

n. t An* = | ni n fe (Aq ik - Aa ik ) - \ (n,An l ) 2 . (A2) 

Next, contracting Eq. (Al) on |cr ifc we get 

\a ik (Aq ik - Aa lk - ArMn fe ) = , (A3) 

and so addition of (A3) to the RHS of (A2) gives 

u t An l = \U lk {Aq lk - Aa lk ) - \h tk Arf An k . (A4) 

Equating the righthand sides of Eqs. (A2,A4), we then have 

n,n fc (Aq* k - A/) = h lk (Aq* k - Aa ik ) + (n,An*) 2 - h tk An*An k . (A5) 

Finally, we contract (Al) on n, and insert Eq. (A2) into the result, thereby obtaining 

An k = rij (Aq ik - Aa ik ) - |n fe n inj (Aq ij - Aa ij ) 

+ \n k (mAr?) 2 - niAn l An k . (A6) 

Let us now turn to the expression (k — k) . Straightaway, we have 

k - k = -DjAn* - n^r'H - An l Ar^ , (A7) 

where the difference in Christoffel symbols is AT l u = T l u — l\. Now, standard formulae 
show both that Al\ = \h kl DiAh k i and h kl D k Ahu = —huD k Aq kl , and with these we rewrite 
the last equation as 

k — k— \n l h kl (D k Ah u - D t Ah kl ) - DjAn* 

+\u l h a D k Aq kl - IhUAn'DiAhu . (A8) 

Focus attention on the next to last term in this equation, which we rewrite as 

^h u D k Aq kl = ln l Ah a D k Aq kl + ^D k Aq kl . (A9) 
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An appeal to Eq. (Al) then shows that 

±n%D fc A g w = \n l Ah a D k Aq kl + ±n;D fc Ao- ifc 

+in,D fc (An l An k ) + n t D k (n< ! A n fc >) . (A10) 

We have seen before that n, and n« are proportional (as both are proportional to dis) and so 
rijA(T lA: = 0. Hence, we may freely shift the derivative D k off of Aa lk and onto rv Next, for 
the last term on the RHS of Eq. (A10), we use r^D^n* = (following from the normalization 
of rf) and find 

^h u D k Aq kl = IrfAhuDkAq 1 * 1 - ±A^ fc D fe n, 

+ ±n,D fc (ArMn fc ) + n fc D, (nW) + D k An k . (All) 

Finally, plugging (All) into (A8), we have the first identity (22a). 

To obtain, the second identity (22b), first simply multiply the first identity (22a) by the 
smearing function N, thereby obtaining a new equation which has Nn k Di(n^ l An k ^) as one of 
its RHS terms. Let us reexpress this term to get the result. Straightaway, 

Nn k D t (n v An k] ) = n fe D* (Nn [l An k] ) + \n k {x\ k An l - n' An k ) D t N , (A12) 

and into this equation we twice substitute (A6), thereby reaching 

Nn k D t (n v An k] ) = n k D t (Nn [l An k] ) + \ [n i (Ag il - Aa il ) - n^n^Ag^ - Aa ik ) 

+n' (m An*) 2 - n t An l An l ] D/JV . (A13) 

Next, appealing to identity (A5), we find 

N^^An^) = n fc D ; (N^Au^) + \ [x\ k \ k Aq a - n% k (Aq lk - Aa lk ) 

+ n'h ifc An i An fc - n;An*An'] D t N . (A14) 

The last equation transforms exactly into 

Nn k D t (n [l An k] ) = n fc D; (Nn [l An k] ) - \ (n k h u Ah tk - n l h lk Ah tk ) D t N 

+\ (n'hifcA^An* - r^Ara'An 1 ) D t N + ^h^Aa^N , (A15) 

and using this result in the new equation obtained from (22a) via multiplication by iV we 
get the desired identity (22b). 
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2. Asymptotic solution to the embedding equations 



We now show how to formally solve the system (32), although our discussion also pertains 
to the system (33) as well as a third system at the next order. As we pointed out in Section 
III.C, the three systems are formally the same. Let us rewrite (32), eliminating / and 
making the substitution g = sin 9g, as 

■ n d 9 dh / A1ft v 

de~d4> = ( A16a ) 

sin !>L + sin 2 e ®L = 7 . (Al6b) 
8<p 89 

Applying 8/89 and (1/ sin 2 9)d/d<p to the first and the second lines of (A16) respectively, 
we find upon elimination of h that 



1 3 / . J\ Id 



sin6>— + 



sin 9 89 V 89 J sin 2 9 8<f> 2 



12 



9 = ii>(B,<f>) (A17) 



with 



iJ ~sin9\89 89J + sm 3 9 8<j)' { } 

Whence we have transformed the system (32) to (A17), that is to say the Poisson equation 

Ag = ^j, (A19) 

where now and from now on A is the Laplacian on the unit sphere. We may similarly 
establish that 

Ah = x , (A20) 

with 

1 fd(3 8a\ 1 8 / 7 \ 
X ~ sin 2 9 \80 ~~84>) + M89 \m) ' ( ' 

Solvability of these S 2 Poisson equations is ensured if the following compatibility conditions 

are satisfied: 

j dQ^j = = j alQx ■ (A22) 

These conditions simply state that the sources ip and x lie i n th e range of the Laplacian. A 
careful treatment of this issue would relate these compatibility conditions to the regularity 
of a, [3, and 7. Indeed, in the axially symmetric case these conditions are tantamount to 
the statement that the B metric o~ a b is free of conical singularities at the north and south 
pole at leading and next-to-leading order. In the general case we expect that the conditions 
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(A22) are also related to the absence of canonical singularities at the poles. Recall that our 
discussion here is also meant to address the system (33) as well as a third system at the 
next order. Such compatibility conditions will also crop up when examining these systems. 
We shall require whatever conditions [on the coefficients A, £>, and Q for the system (33) 
and on similar coefficients for the next order] are necessary in order that these systems are 
solvable. 

A solution to (A19) is formally a solution to the original system (A16). To verify that 
this is indeed the case, integrate Eq. (Al6a) to reach 

h= [ d(f) (sm9% + f3-a J . (A23) 



89 

Substitution of this equation into Eq. (Al6b) then yields 

J dcp sin 3 9 (Ag - ip) = , (A24) 

which clearly holds for solutions to (A19). 

We solve (A19) via standard methods. We expand ip with respect to spherical harmonics: 

^X^WM), (A25) 

l,m 

and are assuming that as a function ip at least lies in the space /^(S* 2 , dQ) of square integrable 
functions. But this is certainly the case, as we expect that the metric coefficients a, j3, and 
7 are C 2 (S 2 ). As discussed below, ip is of even parity and, moreover, must satisfy (A22). 
Therefore, the summation in (A25) is over even values of I and with ^ o — as the constant 
mode so that / is never zero in the sum. 
Search for a solution taking the form 

g = J29i™Yim(0,<f>)- (A26) 

l,m 

Substitution of (A25) and (A26) into (A19) determines the sought-for solution as 

S = -£]7^*WM), (A27) 

l,m 

since, of course, AY lm = —1(1 + l)Y im . Note that in this equation and the ones to follow, 
we never divide by zero, since the constant mode tp 00 vanishes. Similar steps lead to the 
expansion 

h = -Y,j{Tfi) Ylm{dA) - (A28) 

l,m 
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From Eq. (A27) 



g = - sin 9 inTi) Ylm ( d > ® 



1(1+1) 



(A29) 



Whence (32) gives the desired final result, 



f 



a 



Lm 



1(1 + 1) 



d 



cos6Y lm (6A)+^6— Y lm (0A) 



(A30) 



The solution must not to destroy the approximation scheme. We expect a good solution 
provided that the metric coefficients a, (3, and 7 both are C 2 (S 2 ) and lead only to small 
corrections 3 1Zr~ 3 to the leading 2/r 2 behavior for the B Ricci scalar 1Z. 

Let us now consider the action of the parity operator V. Recalling its defining action 
V(x k ) = —x k , one can show that V(r , 9 , <fi) = (r ,n — 9 , <fi + n) and that in (28) the terms 
a and (5 are of the even parity while 7 is of odd parity. Further, we can state that cos#, 
(9/(90, ip from (A18), and A are of even parity, while sin^ and d/d9 are of odd parity. Then 
Eq. (A19) determines V(g) = g, and, consequently, g is of odd parity. In tandem with (32) 
this shows that / and h are of even parity. We now re-express the transformations (30) 
in Cartesian rather than spherical-polar form, adopting standard notation x 1 = x, x 2 = y, 
x 3 = z, and likewise for the X k : 



X = x + r~ 



= y + r~ 



Y 

Z = z + r~ 



xf -yh + 
yf + xh + 



xzg 



\/x 2 + y 2 

yzg 
\Ac 2 + y 2 



+ 0(r- £ ) , 



+ 0(r~ £ ) , 



zf ~ 9\/x 2 + 



+ 0(r~ 



(A31a) 

(A31b) 
(A31c) 



By inspection we see that all middle terms involving square brackets are 0(r°) and of odd 
parity. Therefore, the transformation (A31) may indeed be rewritten in the form X k = 
x k + °C i k + 0(r~ £ ) with an 0(1) odd parity function of the angular variables. Via similar 
analysis we can in principle solve the other two systems we have mentioned [one (33) at 
0(r _1_e ) involving A, £>, and Q, and the other — not written down — at 0(r~ 2 )]. Whence, 
in more detail we find that (34) holds. 
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APPENDIX B: 9 OPERATOR AND FRAME AND CONNECTION 



1. S operator 

Here we develop some standard formulae (drawing mostly from Refs. [29, 30]) necessary 
for and in the notation of this paper. 

a. General two-surfaces 

Suppose that B is a Riemannian two-manifold equipped with metric a a i> and compatible 
covariant derivative operator d a . Let m a , with complex conjugate m a , be a complex null 
vector field on B, chosen such that m a m a = 1. Then the set {e 3 = m, e~ 4 = m} is a complex 
null dyad. (We use 3 and 4 rather than 1 and 2 for the name indices, as we might assume 
our null dyad completes a Newman-Penrose null tetrad. [31] Were we to work with such a 
spacetime tetrad, we would respectively use ej and for the outgoing and ingoing (real) 
null normals to B.) Hatted indices a, b, ■ ■ ■ are null dyad indices, whereas a, b, ■ ■ ■ are general 
frame indices. The complex connection coefficient 

uj = = t\ e% d c e^ a = m a m c d c m a (Bl) 

is the only one which need be considered. Note the the coefficients are not the connection 
coefficients T a b c (which are Christoffel symbols if a, b- ■ ■ are coordinate indices). Introduce 
the operator 5 = m a d a . From Eq. (Bl) we then have 

5m a = ujm a (B2a) 
Sm a = -um a . (B2b) 

Also introduce a dyad derivative operator d a which "sees" dyad indices, and whose action 
on a dyad vector v a is defined as 

d a v £ = e a [v e }+v b u% a , (B3) 

where e a = d a if a, b ■ ■ ■ are coordinate indices. The operator m a d a is known as 3 in the 
compacted spin coefficient formalism. [30] 
Consider a scalar defined on B as follows: 

77 = m ai ■ ■ ■ m a ^m Ci ■ ■ ■ m Cp T ai ... ap+sCl ... Cp , (B4) 
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where T ai ... ap+sCl ... Cp is a rank-(2p + s) tensor field on B (note that s may be negative). Under 
the mapping 

m a ^ e^m a , (B5) 

77 behaves as 

77 i-> e^T? . (B6) 

We say that 77 has spin-weight s, or symbolically sw(n) = s. Notice that 77 is just a 
component 

T s!^si^ ( B7 ) 

p+s p 

of the tensor T ab ... d with respect to the null dyad. More generally, we may define a scalar of 
spin-weight s by taking any rank-g covariant B tensor and contracting any \{q + s) of its 
indices on m a and the remaining \{q — s) indices on m a . For example, T 343 = m a m b m c T a b c 
defines a spin- weight 1 scalar. Of course for our construction here + should be positive 
and an integer. For example, with the tensor T abc one can not obtain a component of spin- 
weight 2 or —2 (components of spin-weight —3, —1, 1, and 3 are possible). Just for the 
sake of concreteness, let us continue the development of the formalism with a tensor T ab ... d 
having the specific index structure given in Eq. (B4). 
Now define the "eth" and "eth-bar" operators, 

877 ee m ai • • • m a ^m Cl ■ ■ ■ m c -m b d b T ai ... ap+sCl ... Cp (B8a) 
877 = m ai • • • m a ^m Cl ■ ■ ■ m^m b d b T ai ... ap+sC1 ... Cp . (B8b) 

By inspection we see that SW(5r7) = s + 1 and sw(5r7) = s — 1. Quick calculations using 
integration by parts followed by appeals to Eqs. (B2) and their complex conjugates show 
that 

<5rj = 5n — swq (B9a) 
877 = 5r] + sujTj . (B9b) 

Moreover, the first equation in (B8), for example, may also be written as 

ST? = m a d a T 3... 3 4 .. 4 . (BIO) 

P+s V 

LEMMA. The commutator of S and 8 is 

(88 - 58)77 = |s^77 , (Bll) 
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where 1Z is the Ricci scalar of B. The lemma is nothing more than the standard "Ricci 
identity" obeyed by covariant derivative operators. In lieu of a proof, we establish the 
identity for a particular illustrative example. With the tensor T abc , build rj = T 343 , a scalar 
of spin-weight 1. Write 5 = m a d a = <i 3 and similarly 5 = m a d a = d 4 . Then we have 

(Bs - 58)?? = y;, ;;i 

= ^334^243 + ^ C 434^~3c3 + ^ C 334^~34c 
= 7Z g^Tg^g + 1Z 434T343 + 1Z 334T343 

= ^"4334^343 + ^3434-^343 + ^4334^~343 ' (^12) 

where TZ^icd are the dyad components of the B Riemann tensor. The dyad curvature TZaidd 
is antisymmetric in its first (and last) pair of indices, and the Ricci scalar is 1Z = — 27^ 34 3 4 . 
Hence, we indeed have 2 [8, = jTZr/. 



b. Round spheres: moving frame and coordinates 

Let us consider a round sphere of radius r sitting in Euclidean three-space E 3 . To 
highlight the fact that we are now working with a round sphere, let us use ji and \x in 
place of m and rh for legs of the complex dyad. On E 3 we choose the moving frame 
{ei-, e 3 , e 4 } = {u, /j, p,}, with v = d/dr and the complex leg 

tfd/dx 1 = -yflr^e-^d/de + i(sin^)- 1 9/90] = ^r^Pd/dC . (B13) 

Here ( = e^cot(9/2) is the stereographic coordinate and P = 1 + CC These conventions 
agree with Dougan's. [32] Clearly, // points everywhere tangent to the foliation of E 3 into 
level-r spheres. Restriction of // to a particular level-r sphere determines a jj, a as before 
(but before denoted m a ). With this choice of fi a we find the following connection coefficient 

^^"vK/r. (B14) 

When B is the unit sphere S 2 , we use the notation lu for this connection coefficient and 
also sometimes a = ^/§C for Here the notation is seemingly odd, but we may adopt 

it to have agreement with the standard Newman-Penrose formalism. [31, 32] Adopting this 
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notation, we may consistently define 

Q r] = (S + 2sa°)r] (B15a) 
Q r] = (S - 2sa°)r] , (B15b) 

where 5o = \j\Pd/d(,. Now the commutator considered in the lemma above is 

(8 c\) - 9o§o)?7 = sr] , (B16) 

since the Ricci scalar for the unit sphere is the constant function 2. Two other useful 
round-sphere formulae are § dil x^oV — ~ $ d£lr}Qox, where sw(x) + sw(rj) = — 1, and 
dkV k = (2V U + (5 Vp + doV^r^ 1 + dV u /dr, where V k is some Cartesian vector. 

c. Spin-0 spherical harmonics 

Let us now document some standard results concerning the unit sphere Laplacian A. 
First, when acting on scalars of zero spin- weight, 

d 



28 5 = t 1 



Sin9 d6 



1 0^ 

+ ^TeW (B17) 



sin eoe 

The eigenfunctions of A are the spherical harmonics Yi m (9,4>) (m runs from —I to / in 
integer steps), which we may also view as functions Yi m ((,0- These obey Sw(Y/ m ) = 
and AYi m = —1(1 + l)Yj m . We may realize the spherical harmonics with the following 
construction. 

Consider the standard stereographic projection. [33] Namely, given a complex number (, 
we produce a unit vector v l = (u 1 , u 2 , z/ 3 ) via the formulae 

u 1 P = ( + C, u 2 P = -i((-(), u 3 P = P-2, (B18) 

whence we indeed have SijV 1 ^ = 1. In these and following formulae, ■ ■ ■ are E 3 Cartesian 
coordinate indices. Viewing ( as a complex coordinate on the unit sphere, we may then 
produce a Cartesian point x % = r(z/ 1 , z/ 2 , z/ 3 ) given "polar coordinates" (r, (). The Cartesian 
components, 



// = y/Ir-ipdx'/dt = ^J\P~ l (l - C 2 , -i(l + C 2 ),2C) , (B19) 

of the vector (B13) are obtained via the chain rule. Now choose an arbitrary constant 
Cartesian vector field c\ Clearly then sw(//q) = 1 and sw(z/q) = 0. It follows that, 
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viewed as functions on S 2 , the components ji l and v % have spin-weights 1 and respectively. 
Therefore, computations based on the rules (B15) give the useful formulae 

g z/ = , g o/ / = , %tf = -z/ . (B20) 

With these formulae, it is easy to verify that 26 9o^ = — 2i/; that is to say, the v l are 
essentially the Y lm . Moreover, defining f lj = v % v i + + fi 1 ^ : , one can check with (B20) 
that 

2oWV - IP) = -6(i/V - IP) , (B21) 

thereby showing the z/V- 7 — |/ u to be essentially the Y 2m . Naively, there are six z/V J — 
but only five Y 2m . However, notice that v l v i — = (one condition); hence, z/V J — |/ u has 
only five independent components as expected. As is well-known, one may continue building 
higher zero spin-weight harmonics as symmetric trace-free Cartesian tensors polynomial in 
v % and f lJ . 

d. Spins spherical harmonics 

Extension of the action of A = 2§ 9 to scalars of arbitrary spin-weight is trivial, since 
25 9o is defined on arbitrary scalars. Eq. (B16) shows that on an s spin-weight scalar 

25 Bo = A - 2s . (B22) 

The spin-s spherical harmonics s Yi m are the s spin-weight eigenfunctions of the Laplacian 
A. Following Goldberg et. al. in Ref. [34], we define 

s Y lm = [2 s (l - s)\/(l + s)!] 1/2 5^ m (B23a) 
- s Y lm ee [T(l - s)\/(l + s)!] 1 /2(_i)-g-y Jm , ( B 23b) 

here with the restriction < s < I. The discrepancies in factors of 2 s / 2 with the definitions 
in Ref. [34] stem from the fact that the 5 operator of Goldberg et. al. is y/2 times our 
own. Notice that we are not allowed to increment the spin-weight beyond the range —I to 
/ (the same restriction on the integer m). This is suggested by the considerations above. 
Indeed, the formulae in (B20) clearly show that both 5q and Sq annihilate is 1 ; that is to say, 
Soi^im = = So-i^im- Furthermore, repeated use of these formulae shows that both 5q 



44 



and 5q annihilate z/V- 7 — our geometric representation of the Y 2m . This then implies 
So2^2m = = 8 -2^2m- These concrete examples exhibit the idea behind the identities 

9oi^m = = 9 -j^. (B24) 

The standard identity Y/ m = (— l) m Yi m along with the definitions of 9o and 5o determines 
that 

s Y lm ={-l) m+a _ s Y l> _ m , (B25) 

now with —l<s<l. We have the main 

Lemma. With the Laplacian A = 23 3 and for —l<s<l, 

A s Y lm = [s(s + 1) - 1(1 + l)\ s Y lm . (B26) 

To prove the lemma, first establish the positive s case via a simple induction argument, one 
using Eq. (B16). Next, obtain the negative s case with the identities (B25) and (B22). 
With the formulae collected so far, we gather the results, 

3o,Y lm = [(l-s)(l + s + l)/2]^ 2 a+1 Y lm (B27a) 
Qo.Y lm = -[(I + s)(l-s+ l)/^ 1 / 2 ^^ , (B27b) 

which augment those given in Eq. (B24). Hence we may view So as a kind of raising operator 
and 5 as the corresponding lowering operator. 



2. Frame and connection 



The spatial frame 





1 (d 


dx a J 


" M [ik 


<• 9 


1 


dx a 


\m,Q 2 — 



e h = -fl ( 4- - W a ^-] = ±-(^- Wrnm a ^- - W m m a ^- ) (B28a) 
M\dr dx a ) M\dr dx a dx a ' 

d d 

8(- mc SC 



is dual to (63). Note that fn^ = rh^ and = fh^. The basic coframe variables are then 
seen to be M, W m , m^, and m^. 

We may now use the method of Cartan to calculate the connection coefficients de- 
termined by the co-frame (63). The method starts with the no-torsion formula 

de^ = -J 1 1 A e l , (B29) 
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and the identities 
dr A d( 



1 



M(|m c -| 2 - |m f | 2 ) 
1 



m^e h A e 3 — m^e h A e 4 



) (B30a) 

dC A dC = M(|m . | ;_ K|2) (Me 3 A e 4 + A e 3 - A e 4 ) (B30b) 

prove useful in carrying out the necessary calculations. In any case, the list of u>0 is the 
following: [cf. Eq. (86)] 

Wa _ h = -SlogM (B31a) 



^433 



W 3h4 



\m? 



|m c |< 



1 


"(1 




2 _ 




2 y 


2M 






2 _ 


m c 2 



3h3 



w 34h - 



M(\m- C \ 2 - |m c | 2 ) 



(9m^ drriQ 



(B31b) 
(B31c) 
(B31d) 



2M(|m ( -| 2 - |m c | 2 ) 



1 



+ M ~ + ^^433 - W m UJ m ) , 



(B31e) 



where the prime denotes differentiation by d/dr. Moreover, with u = w 43 3, the action of 9 is 
that of e 3 , e 3 — a;, eg + u, on log M, W m , and W 7 ™ respectively (consistent with the defining 
action of 5 given above). 

Let us obtain an expression for the scalar curvature 1Z of B. Denote by 9 3 = fn a dx a the 
pullback of the co-frame e 3 to a B surface, where d is the B exterior derivative (the exterior 
derivative of X has appeared as d). Then the connection one-form on B is 



^433# S + ^lO l = W0 S - UQ 1 . 



'434 L 



(B32) 



In terms of the dyad components of the B Riemann tensor, the B scalar curvature is 1Z 
—27^4343, whence we may obtain 1Z from the formula 



^43430 A ^ 



Im 



d(u6*) 



We find 



Tl = -2(e 4 [cj] + e 3 [cu] + 2uu>) 
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(B33) 



(B34) 



as the result. 

APPENDIX C: DIVERGENCE OF THE BASE HAMILTONIAN 

Assuming both the fall-off given in Eqs. (7,8) and the center-of-mass scenario with 
lapse given by Eq. (11), in this appendix we examine the "base Hamiltonian" ifs, isolating 
its divergent contribution in the r — > oo limit. We point out that our calculations here 
complement those given in Appendix C. of Ref. [4], which examined the same issue via a 
different method. The method we adopt here brings our Section IV results into sharper 
focus. Moreover, the discussion in Section Vis based in part on the results of this appendix. 
We now consider if s off-shell; we do not assume that the £ initial data obeys the scalar 
constraint H. — 0. 

1. Geometric identities and two lemmas 

Before turning to i?s, let us first collect some geometric identities and prove two lemmas. 
Consider the following expression for the E scalar curvature: 

R = lZ J r k 2 - k ab k ab + 2Dj (kn 1 + V) . (CI) 

In this equation we may write k 2 — k a bk ab = \k 2 — 2k mm kmm with k = 2/c mj7i as before. 
Moreover, here bj = n k D k nj = —a k D k \ogM, where the last equality follows from the fact 
that the £ Levi-Civita connection is torsion-free. One way to verify it is to use rij = MDjr 
in n k DkUj and then carry through with some algebra, along the way using the identities 
DkDjr = DjD k r (by no torsion) and n k DjUk = (by the normalization of n k ). Note that 
V points everywhere tangent to the B foliation (i.e. bjn^ = 0), whence the single complex 
component b m = — 31ogM completely specifies this B vector. A simple strategy based on 
the Ricci identity for deriving (CI) is outlined in the appendix of [5]. From Eq. (CI) we get 

R = <Jl-\k 2 - 2k mm krn m + 2M~ 1 k' - (W m Bk + Wrffik) - 4M _1 S§M . (C2) 

To reach this equation we have used the following identities: Djn? = —k; DjV = o^DJoj + 
n l n j Dibj = (db^ + db m ) — b a b a , where the term within the parenthesis is the B divergence 
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of the B vector b a ; and 

n k d/dx k = M~ l (d/dr - W m B - W^) , (C3) 

which is just Eq. (B28a). In Eq. (C2) the terms involving k mm kmm and W m dk are each 
0(r~ 4 ) and of even parity, as shown by the asymptotic expansions (67b, 76, 77,78) and the 
parity discussion found in Section IV. B. Moreover, from Eq. (36c) we infer that 

M" 1 ~ 1 - x Mr- x - (1+£) Mr" 1_e - ( 2 M + E.P.T.)r~ 2 , (C4) 

where we recall that l M is of even parity. These considerations and our Section IV results 
show that R has the expansion R ~ 3 Rr~ 3 + ( 3+£ ) Rr~ 3 ~ £ + 4 Rr~ 4 with 

3 R = 2 2 k + 3 IZ- 4(1 + do^M (C5a) 
^R = 2(1 - e) {2+£ ^k + ^IZ - 4(1 + 5 5 ) (1+£) M (C5b) 
4 i? = ±n - 4(1 + 5 B ) 2 M + E.P.T . (C5c) 

The coefficient 3 R is easily seen to be of even parity via Section IV results and arguments. 

We now state the lemmas. Let fl^ v = fi^iV 1 as before. Then we have the following: 
Lemma 1. 

j dn/3\ (3+£ ^R ={l-e)j d£l/3\ (3+£) [o-^a x ^ XuK ] . (C6) 

Lemma 2. 

fdnp ± v A R = 0. (C7) 
To prove Lemma 1, first use Eqs. (96,100a,C5b) to obtain 

jdnp\^R = {l-e)jdnp\^[o ii o* l R ikjl ] 

+ j dVlp ± u [e {3+£) IZ - 4(1 + g 8 ) {1+£) M] , (C8) 

and then with Eq. (98) replace the £ Riemann tensor on the RHS with the spacetime Riemann 
tensor. The second integral on the RHS vanishes. Indeed, double use of angular integration 
by parts along with Eq. (B20) shows that the unit-sphere average of P ± v (l + <3o3o)^ 1+e ^W 
is zero. Moreover, the unit-sphere average of (3 1 - U ( 3+£ ^1Z is also zero. Indeed, explicitly 

^IZ = + B 2 ^ - 23 g 6^ - 2b, p , (C9) 
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that is to say the result (79) only with a^- replaced by This result can be verified with the 
techniques of Section IV. Angular integration by parts along with (B20) again shows that 
the average in question vanishes. (Note, however, that the unit-sphere average of ( 3+£ ^lZ 
need not vanish. The discussion in Section V of the mass-aspect for the energy scenario 
rests on the fact that the same is true for 3 1Z.) As for Lemma 2, the same arguments apply, 
because we expect 4 1Z to have the form 

A U = & Q c m + 8g Cwt - 28 5oc M/i - 2c^ + E.P.T. (CIO) 

This result should follow from Eq. (C9) upon setting 5=1 and including quadratic even- 
parity corrections. 



2. Divergence of the base Hamiltonian 

Let us now recall that H% is the volume integral given in Eq. (6), where 

n =h( Ki > Kij - K (C11) 

In what follows, let us assume for simplicity that the slice £ is topologically Cartesian three- 
space R 3 , that we may ignore the issue of inner boundary terms. Upon integration against 
the product of the £ volume element d 3 x\f~h = d 2 xdr^faM and lapse N, the terms in (Cll) 
which are quadratic in the £ extrinsic curvature tensor do not contribute to the limit. 
Indeed, these terms are 0(r -4 ) and of leading even parity, the volume element is 0(r 2 ) and 
of leading even parity, while iV is 0(r) and of course of leading odd parity. Then the product 
of all of these terms is 0(r~ l ) and upon integration over the radial coordinate would yield 
a logarithmic divergence were it not for the fact that this product is of leading odd parity. 
Let us then focus on the term 

1 f .„ rr I 

16tt 

where we have used Eqs. (ll,36c,40) and the parity properties of leading terms to isolate on 
the RHS only those terms from the LHS which lead or could lead to an infinite limit. Next, 
we insert the radial expansion for R and integrate in r term-by-term, thereby reaching 



/ d 3 xv / hNR — [ dmrr 3 (3 L v R, (C12) 

JS 167T Jt, 



16tt 



16tt 



d 3 xVhNR ~ (C13) 

j dni3 L u 3 R - m ^ £ _ - j dnp\^R dnp\ A R . 
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The first term on the RHS vanishes due to the even parity of 3 R. These calculations, the 
argument of the previous paragraph, and the lemmas above then show that 

Hx ~ j <mp\ (3+£) [a^a x ^ XuK ] + • logr . (C14) 

Note that the coefficient of logr has been found to vanish both by parity arguments and 
Lemma 2. Because the first divergent term in Eq. (C14) has the opposite sign from the 
corresponding term in the main text [cf. Eqs. (5,11,60)], we see that the full Hamiltonian 
H-z + Hb is finite even off-shell. Whence Hb must be finite on-shell. 
Finally, as an aside we note that the results of this appendix show 

M^ k = I dtl 4 [nViy v k (C15) 

to be equivalent to Eq. (16) as a definition of center-of-mass coordinates for initial data 
sets. This expression in tandem with the GauB-Codazzi-Mainardi splitting of the spacetime 
Ricci tensor might be used to derive other equivalent expressions for center of mass. 
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